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Abstract 

Given that observations seem to favour a £!o < 1, corresponding to an open 
universe, we consider gauge - invariant perturbations of non-flat Robertson - Walker 
universes filled with a general imperfect fluid which can also be taken to represent 
a scalar field. Our aim is to set up the equations that govern the evolution of the 
density perturbations A so that it can be determined through a first order differential 
equation with a quantity /C which is conserved at any length scale, even in non-flat 
universe models, acting as a source term. The quantity JC generalizes other variables 
that are conserved in specific cases (for example at large scales in a flat universe) 
and is useful to connect different epochs in the evolution of density perturbations 
via a transfer function. We show that the problem of finding a conserved K. can be 
reduced to determining two auxiliary variables X and Y, and illustrate the method 
with two simple examples. 
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1 Introduction 



Despite much theoretical progress made in the last fifteen years, the problem of large- 
scale structure formation in the universe remains an unresolved puzzle for cosmologists. 
In fact, as new and more reliable data becomes available, none of the models proposed 
during this time seem satisfactory. 

In the standard Hot Big Bang model jjj , the most popular mechanism advocated to 
explain structure formation is that of the gravitational instability: small perturbations 
created during an inflationary era eventually grow in the matter dominated era. It is 
therefore necessary, in order to use observational data to constrain theoretical models, 
to relate the original amplitude of a perturbation mode at an initial time to that of the 
same mode today, i.e. to derive a transfer function (cf. @, §)• 

In cosmological perturbation theory the evolution of a perturbation mode is usually 
given by a second order ordinary differential equation (e.g. see [|], || [| |2| ) . It is therefore 
very useful to set up the problem in such a way that, at least in the simplest cases, one 
has a conserved quantity at ones disposal, so that the problem is effectively reduced to 
solving a first order differential equation with a constant source term which, at least in 
principle, can always be integrated. 

Indeed, such conserved quantities can be found in flat (see @, |, | and g(], 0, g|, |; 



and almost flat [13, [Tl|]) universe model, for a perturbation mode whose wavelength is 
much larger than the Hubble radius. This is quite satisfactory in the standard inflationary 
scenario where the present day density parameter Oo is predicted to be equal to 
unity with a high degree of accuracy. However, there is no compelling observational 
evidence for the case of Oo = 1; rather, a low density universe seems to be favoured by 
analysis of present day observational data [15]. 

For this reason, and taking into account the appeal of the inflationary paradigms and 
the fact that this is up to now the only scenario in which perturbations are predicted 
(by quantum fluctuations within the horizon) rather than assumed as initial data, some 
authors [16, 17, 18, 19 1 have considered the possibility of an era of "minimal" inflation, 
in which the accelerated expansion does not last long enough to drive Qq to unity today. 
This possibility can be a source of a certain amount of debate [2C], raising on one hand 
arguments of fine tuning (which however is in one way or an another needed for other 
reasons in most of the inflationary models proposed up to now), and on the other hand, 
"how to put a measure on the set of initial data" (lacking this measure, according to this 
point of view, the fine tuning argument is a weak one). 

In this paper we do not enter in such a debate, instead we take the case of a negatively 
curved universe as an open possibility, and examine perturbations in this context. For 
the sake of generality, we shall consider an imperfect fluid, which later we shall take to 
represent a perturbed minimally coupled scalar field. 

Our aim is to set up the equations that govern the evolution of the density perturba- 
tion mode A^ so that it can be determined by a first order equation with a conserved 
quantity lO- 1 ^ as a source term. In other words, we define a simple algorithm that in 
principle allows us to find a new quantity KS 1 ^ which is conserved at any scale, even 
in non-flat models. Then K, generalizes known conserved quantities, and allows one to 
define a transfer function in a more general set of circumstances. 
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As it will be shown, the problem of finding the quantity KS V ' in the cases of: i) 
adiabatic perturbations in a perfect fluid or ii) perturbations in a minimally coupled 
scalar field, reduces to determine two auxiliary functions X and Y. 

It is noticeable that a completely consistent theory of large scale Cosmic Microwave 
Background (CMB) fluctuations in open models seems to be lacking until now 21]. 
With the cautionary note that even the concept of power spectrum has not yet been 
generalized to include perturbations over super horizon scales in open universe models, 
we hope that the work presented here is a step towards a consistent and useful treatment 
of fluctuations in f2o < 1 universes. 

This paper is organized as follows: In section 2 we discuss very briefly the formalism 
and variables used. In section 2 and 3 we define two basic curvature variables and discuss 
under what conditions they are conserved both in the case of an imperfect fluid and a 
minimally coupled scalar field. In section 4 we give a prescription for finding quantities 
that are conserved at all scales for any given background evolution. In section 5 we apply 
our algorithm in two simple cases: i) a coasting (f2 = const) phase and ii) a De-Sitter 
inflationary phase. In section 6 we briefly summarize how to use our algorithm in order 
to find a transfer function. Finally in section 7 we end with a discussion. 

In this paper we assume standard general relativity, with Einstein equations in the 
form R ab — T;Rg a b + h~9ab = K T ab , K an d A. being the gravitational and cosmological 
constants respectively. 



2 Formalism 

For the sake of self- consistency, we sketch in this section the covariant approach to den- 
sity inhomogeneities in a general curved spacetime [p2|] . This as been applied to almost 
Friedmann - Lemaitre - Robertson -Walker (FLRW) universes dominated by a barotropic 
fluid |23|, 0], as well as by a simple scalar field a general imperfect fluid pO|, [| 

and a mixture of interacting fluids [24, 25]. Here we just give an outline of the material 
which will be needed later on, with emphasis on known conserved quantities. 

2.1 Covariant cosmology 

Let us consider an imperfect fluid flow in a general curved space -time. The 4 -velocity 
of the fluid is the tangent along the fluid flow: u a = dx a /dr (u a u a = — 1), where t is the 
proper time along the flow lines (the world lines of observers comoving with the fluid). 
Introducing a projection tensor orthogonal to [26] h ab = g a b + u a Ub (h a b Ub = 0) the 
energy - momentum tensor is decomposed as follows: 

T ab = HU a u b + ph ab + 2q {a u b) + ir ab , (1) 

where \i and p are the energy density and pressure of the fluid, and q a and 7T a b are the 
energy flux and anisotropic pressure contributions. When the fluid is not perfect (i.e. in 
it's perturbed state) the fluid 4- velocity u a is not uniquely defined M, |25| : it is however 
standard to choose u a either as the particle frame, in which T a b has the general form (|l|) , 
or as the energy frame, in which q a = (for a discussion of these choices, see |]]). Here 
we shall take the particle frame choice, as this allows one also to describe more general 
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cases for example a non - minimally coupled scalar field and higher order gravity theories 

13- 

The covariant derivative of u a can be split into four parts: 

u a -b = \®h ab + a ab +u ab - a a u b , (2) 

where G = u a - a is volume the expansion, a ab = ai ab \ is the shear (a ab u b = a a a = 0), 
oj ab = mum is the vorticity (u> ab u b = 0) and a a = u a = u a - b u b is the acceleration (the dot 
denotes the proper time derivative). It is useful to introduce a length scale factor along 
flow lines by the relation 

when the universe is an exact FLRW space -time H is just the usual Hubble parameter. 
The evolution equation for the expansion is the Raychaudhuri equation 

9 + ±6 2 + 2{a 2 - to 2 ) - a a - a + + 3p) - A = , (4) 

where a 2 = \o ab o ah and to 2 = ^u> ab u> ab are the shear and vorticity magnitudes. 

The matter equations of motion T ab . b = are equivalent to the energy conservation 
and momentum conservation equations: retaining only first order contributions in the 
inhomogeneity variables (the reason for this will be clarified in the next section) we have: 

fi + 3hH + h {3) X7 a V a = , (5) 

and 

ha a + Y a + h [F a + n o ] = , (6) 
where for later convenience we define 

h = (ji+p), *a = j, ,F a = ^ a -(3c 2 s -l)H^ a , U a = ^V b TT ab . (7) 

An exact FLRW can be covariantly characterized by the vanishing of the shear and 
the vorticity of u a and by the vanishing of the spatial gradients (i.e. orthogonal to u a ) 
of any scalar /: 

a = u = , ( 3 )V a / = ; (8) 

in particular the gradients of energy density, pressure and expansion 

X a = ^V af i, Y a = ^V aP Z a = ^X7 a G, (9) 

vanish, where Y a = a a = 0. Then fi = fj,(t), p = p(t) and G = @(t) = 3H(t) 
depend only on the cosmic time t defined (up to a constant) by the FLRW fluid flow 
vector through u a = —t^ a . The energy momentum tensor ([]]) necessarily reduces to the 
perfect fluid form T ab = [iu a u b + ph ab , and it follows that these models are completely 
determined by an equation of state p = p(fJ-), the energy conservation equation (a) (with 
no imperfect fluid contributions) and the Friedmann equations 

3H + 3H 2 + \k{p + 3p) - A = , (10) 

and 

F 2 + f = ±^ + ±A, (11) 

where the latter is a first integral of the former (when H ^ 0), which is (Q) specialized 
to a FLRW model. 



3 



2.2 Fluid inhomogeneities and gauge invariant variables 



The vector Y a is one of a set of variables we can define [22] to describe the fluid inhomo- 
geneity: 

Y a = h a b P;b , V a = £^, Z a = £h a b @. b] (12) 

these are the covariantly defined spatial (i.e. orthogonal to u a ) gradients of the pressure, 
energy density and expansion; more precisely the dimensionless comoving fractional spa- 
tial gradient of the energy density D a is the physically relevant variable to characterize 
the energy density spatial variation. All of these covariantly defined exact variables can 
be in principle measured as spatial gradients by observers at rest with respect to the 
fluid, and in general (i.e. in an arbitrary spacetime) they obey exact non- linear and 
covariant evolution equations which have been derived in for a perfect fluid and in 



[24 1 for an imperfect fluid; these equations are coupled with the exact evolution equations 
for a ab and uj ab (see e.g. |26|). 

The above introduced covariant gradients Z a and D a , and in general any variable such 
as those in (||) that vanish in the FLRW background, are gauge - invariant (GI) variables 
(for more details, see [||, || and |2Sj|), so we can proceed to obtain equations for 
them valid in an almost FLRW universe (a space - time where these variables are small) , 
approaching this universe from a general space -time rather than perturbing an exact 
FLRW model. Then the linearization procedure we apply consists in dropping terms 
such as a 2 in (Q), i.e. terms which are of second order in the GI variables, retaining only 
terms linear in T> a and Z a and using zero -order expressions for fi, p and O = 3H [this 
clarifies the meaning of neglecting higher order terms in (||) and (||)]. 

At zero - order these variables are functions of the cosmic time only, and satisfy the 
usual equations valid in an exact FLRW universe. Applying this linearization procedure 
to a universe where the fluid is barotropic [p = p{fJ>)], one obtains a closed pair of first - 
order (in proper time derivatives) equations that couple V a and Z a only [[2^]; this is 
equivalent to a second - order equation for V a [ 23 , |{| . It was noted in || that in general 



a source term due to vorticity is present in this equation; it also admits an entropy 



perturbation source term when a more general non barotropic fluid is considered 1C ] . 
Finally, we note that our covariant exact variables can be related to those of Bardeen 
at first order (see || |l0|). Using this it is then possible to connect our GI perturbation 
variables to the standard analysis of the spectrum of perturbations and fluctuations in 
the cosmic microwave background, which are usually given in some specific gauge. 

2.3 Conserved quantities 

In order to proceed towards the identification of GI conserved quantities, we introduce 
the 3 -curvature scalar of the projected metric h ab orthogonal to u a ; which in general is 
(see the appendix in ]||) 

^R = 2{-\Q 2 + a 2 -uj 2 + k^ + K) , (13) 

and when u = it reduces to the 3-curvature of the surfaces orthogonal to u a . From 
(13) we can define the GI curvature gradient as 

C a = £ 3 h a b ^R b = -\Qt 2 Z a + 2^£ 2 V a , (14) 
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where the second equality holds in the linear approximation; then this expression shows 
that only two of the three variables T> a , Z a and C a are linearly independent at first order. 

Up to now we have considered vectorial GI variables. However we are interested 
here in characterizing matter clumping through scalar density perturbations: instead 
of introducing these through the standard non-local (ADM) splitting |2g|] , we simply 
locally define GI scalar variables by taking the divergence of the vectorial quantities. In 
particular we have: 

A = lh ah V a . b , Z = £h ab Z a , b , C = -\®£ 2 Z + 2k^i£ 2 A , (15) 

where A is the locally defined density perturbation scalar variable describing matter 
clumping (the factor £ in these definitions comes in for dimensional reasons). We also 
introduce the related quantity 

C = C -4if7 -1 A , 7 = l + u>, w = -, (16) 

which reduces to the GI curvature perturbation C when K = [K is the curvature 
constant of FLRW models appearing in (11)]. Once a harmonic analysis is carried out, 



A is equivalent to the Bardeen variable e rn ]10|, g] and C corresponds to the variable £ 
defined in [p]] (see |i~i"H ). Again only two of our variables are independent, therefore the 
A evolution can be obtained from a linear system of two equations for A and C or A 
and Z. For a general fluid with speed of sound c 2 = dp/dfi we obtain ||: 

A-{3Hw-[f -§} H-^A-^C 
= 3£(l + w)H[F + U] -e(l + w)V 2 V , (17) 



Z + 2HZ + ±K/iA + 



((3) V 



2 | 3K 



-i i^v 2 + m) [f + n] + |/t [f + n] 



(l+w) { V 



(18) 



and 



C 



it 2 Hc 2 s (3)y2A _l M 2 Hw \{3)tj2 , 3K 
(l+w) (l+w) [ v "•" ~W 



+4£ 3 i? (3) V 2 [F + n] + 



AK£ - 2th 



( 3 )V 2 \I' 



(19) 



where 



n = £ (3) V a n a , F = £ (3) V a F a , * = ^ (3) V a * a 



The system of equations (OM) admits 7 ^ (when 7 = one has c 2 



w 



(20) 
— 1) and 

is derived for general K and A; however in the following we shall assume A = 0. 

It is useful at this point to carry out a harmonic expansion of our variables introducing 
scalar harmonics which satisfy the Helmholtz equation ^SJ 2 Q^ = —k 2 /£ 2 Q^ k \ 
where A; is a comoving (i.e. constant) eigenvalue. If v is a non- negative real wavenumber 
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then in a flat K = universe, it is associated with the physical wavelengths A = 2irl/iJ^ 
since when K = 0, k = u, however if K = — 1 then k 2 = v 2 + 1 [29|. Since observations 
tend to favor Qq < 1 115[| , we shall now restrict attention to the cases K = —1,0. We 
consider the equations for A and C, since it is the latter variable (and not the GI 
curvature perturbation C) that will be eventually conserved when Q ^ 1. Writing 
equation (19) in terms of harmonic components, the evolution equation for becomes: 





C 



M 2 H 3 c 2 ( „2 
(l+w) 



K 
WW 1 



AM + 



4£ 2 H 3 w 
(1+w) 



AK 
WW 



v2 
WW 



K 
WW 



- # 2 



4i« - 2r/i 



WW 



(21) 



For large scales, i.e. wavelengths A 3> -ff 1 =>■ <C 1 this equation reduces to 

- W 4£ 2 H 3 c 2 ( K 



° - + (!+„,) ^F^J A + IT 



+w) 



4A - 

FTP' 



+ 4rJT 



3 z_r3 / X 



+ # 2 



4iO - 2£ 3 /i 



(22) 



Let us now examine the consequences of this equation for a number of different cases. 

If the background is flat, K = 0, we immediately see that C = C is conserved in 
the large scale limit even when entropy perturbations and imperfect fluid source terms 
are present. This is an important result, since we can now use this conserved quantity 
to write down a general solution of the first order equation for A^ u \ (0), in the long 
wavelength limit: 



am 



t r rt 



/o 



i A ^dU 



B^(h)dti, 



where 



K,h 
2H> 



B {u \t) 



(i- 



4f 2 H 



A(t) = 3Hw 



1&W + 3i{l + w)H [fM + nM 



(23) 

(24) 
(25) 



and to is the epoch at which an initial condition for A^ is specified. 

The conservation of C is also very important in the study of perturbations in both 
standard and non-standard inflationary models (based on generalized gravity theories^) 
since it can be used to directly connect the amplitude of present day large scale structure, 

*A completely consistent theory of CMB fluctuations in open models seems to be lacking until now 
pl| ], and even the issue of completeness of the set of eigenfunctions of the Laplacian operator seems 
unresolved. 

2 A brief summary of the needed harmonic functions is given in Appendix B of Q. 

3 The fluid flow approach to perturbations can also be used to study perturbations in generalized 
gravity theories. The basic idea is to treat all additional contributions to the field equations except the 
Einstein tensor part as contributions to an effective energy momentum tensor. Effective fluid quantities 
are then easy to compute and GI quantities based on spatial gradients can be defined as usual (see H). 
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which came inside the horizon during the matter dominated era, to the initial conditions 
just after horizon - crossing during the inflationary era. 

We now turn to the case when the background is an open (K = —1) model. In this 
case C is conserved if the following inequalities hold: 

< 1 and -2—2 < 1, (26) 



£ 2 H 2 ^ " £ 2 R 2 

i.e. if the long wavelength limit is valid up to v = 1, provided that the entropy pertur- 
bations and some of the imperfect fluid terms do not exceed A^. However, in an open 
universe the cosmological density parameter f2 is given by 

so the second inequality only holds if O is close enough to 1, but then C is also approxi- 
mately conserved. 



3 Scalar fields 

Let us now extend the analysis carried out up to now to a universe dominated by a scalar 
field (j), as is the case during an inflationary phase; here we give only a brief outline of 
the method (details can be found in fll^l ). We assume for simplicity that <fi is minimally 



coupled (the non - minimally coupled case will be considered elsewhere [31[): then the 
Lagrangian density of <j> is 

^ = -V^g[^' a + v(<p)} , (28) 

where V((f>) is a general unspecified potential assumed to be responsible for inflation. We 
shall also assume that (f) a is a time - like vector [^?J , and we make the choice 

u a = -0- V" . (29) 

The former assumption will be well justified by a continuity argument in an almost FLRW 
universe, as it clearly holds in an exact FLRW model; the latter choice of u a is uniquely 
determined by 4" a on requiring u a u a = — 1 (note that uj a b = with this choice). Then 
the energy - momentum tensor of 4> takes the perfect - fluid form (|l]), with p = ^tp 2 — V(<p) 
and a = \ ft + V (</>). With our choice of u a the Klein - Gordon equation D0 — V'((p) = 
for (j> (the prime indicates the derivative with respect to <j>) takes the form 

4> + + V'(<j>) = (30) 

in any space -time, while the spatial gradient of <f> identically vanishes: h a (j)^ = 0. A 
gauge - invariant variable associated with is h a b (j)^', from this we can define a dimen- 
sionless gauge - invariant scalar variable $ simply related to the density perturbation 
A: 

$ = ih oc {l4r x h a % h ). c A = 7 $. (31) 

It is then easy to show that equation (|l7]) for A still holds, whereas the equation for C 
contains a further term which can be interpreted as due to an entropy perturbation |T^. f?j. 
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The equation for the harmonic component C^"' of C is then 

AM , (32) 



O — 

7 



v 2 6|K| / 




while ( |l~7| ) still holds for the components A^ and C^ u \ Note that we left an explicit 
| if | factor in the above equation to show where K appears: when K = 0, C ([!(]) reduces 
to C and the last term in ( |32|) containing \K\ can be dropped, i.e. scalar field and 
barotropic fluid perturbations satisfy formally similar equations [(^) reduces to ( |l9| ) on 
putting c 2 = 1; however this is only a formal identification: c 2 = c 2 (t) ^ 1 for a scalar 
field]. When K = — 1 instead, the second term in the square brackets in ( |32[ ) will not 
be negligible in general, even for scales A ^> H : C will then decrease if V > 0, but it 
will increase in the opposite case. We can ask now in which cases will be conserved 
for scales far outside the effective horizon if -1 when K = — 1. It is clear from the above 
equation that in general C will not be conserved, even if v 2 /(iH) 2 « 1, as the other 
term in the parenthesis will be not negligible in general. However let us consider the 
horizon crossing condition 

u/(£H) m 1 & v w (1 - > l , (33) 

where Qc is the value of f2 when the comoving scale fixed by v crosses the horizon; we 
see from the above relations that scales with v < 1 never cross the horizon. Scales which 
are at present inside the horizon, i.e. A < Hq 1 , correspond to v > vq ~ (1 — f2o) _5 
and therefore left the effective horizon when 1 — 0,c < 1 — ^o- in standard inflationary 
models that predict |1 — fio| to be exponentially close to zero all scales of major interest 
left the horizon at the very end of inflation, and the corresponding wavenumbers v are 
exponentially large; if instead |1 — Oo| ~ 1 (i.e., if we have a minimal inflation), than we 
can see many of the scales that left the horizon from the very beginning of the inflationary 
phase, and in this case v > 1. In any case these scales, and those outside the present 
horizon (A > Hq 1 ) that lie in the interval 1 < v < (1 — f2o)~ 5 (practically vanishing 
for minimal inflationary models), were far outside the horizon only when Q was already 
close enough to unity, 



A » H- L ^ v{\ - «; 1 ^ I - 0« 1 - n c , (34) 

in the last stage of inflation. It follows that, in a universe which is fairly open today, 
cannot be conserved for these scales in epochs when |1 — Q\ ~ 1. For these scales 
C^ v ' will be practically conserved whenever f2 is very close to unity, provided also 

a, 7 

+ t~1, (35) 



3Hcf) 6 

e.g. during a slow rolling phase implying <j) <C 3H(p. 

Finally, those scales larger than the present horizon Hq 1 corresponding to v < 1 have 
always been outside the horizon, so the condition v 2 /(IH) 2 <C 1 can be satisfied even 
when |1 — ft\ ~ 1; however at these scales CM will be not conserved in general, unless 

1 V 

77 r < 1 . (36) 

6 3H(j) 
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One is mainly concerned with the evolution of cosmological perturbations which are at 
most of the size of the present Hubble radius or slightly larger: if |1 — Oo| ~ 1 pertur- 
bations on such scales cannot be far larger than the effective horizon from the moment 
they exit the horizon itself up to when f2 is close enough to unity. During such an epoch 
a quantity such as C cannot be conserved on these scales, and the density perturbation 
evolution equation (|l7| ) is no longer an independent equation with a constant source term 
responsible for the growth of the density perturbation itself. In the next section we will 
construct general conserved quantities valid for particular inflationary models. 



4 Conserved quantities 



In the previous section we discussed the evolution of the scalar curvature variables C 
and C and showed that for a scalar field neither of them are conserved in general when 
K ^ for scales much larger than the Hubble radius. It is interesting to ask whether 
it is possible to construct a quantity K, from the set of perturbation variables that is 
conserved at any scale, even for It would therefore generalize the two curvature 

variables already discussed. 

Let's start by defining a general scalar quantity that generalizes C and C: we found 
convenient to choose that it take the form 



/C = -±£ 2 ®Z + 2^£ 2 (l 



A + eXZ + eYA, 



(37) 



where X = X(t) and Y = Y(t) are two functions to be determined, and e = 0, 1. For 
e = K, reduces to C, and it reduces to C if in addition we choose K = 0. As we are 
now going to show, for e = 1, there are cases (clarified below) in which we can determine 
X and Y so that K, is conserved at any scale (i.e. any mode K^ u ' is conserved), both for 
a flat or for an open universe (the following equations hold only for K = 0, — 1). To find 
an evolution equation for K, , we take its time derivative and use the equations for Z 
(18) and A (|l7|). Writing in terms of harmonic components, and keeping the imperfect 
fluid source terms for generality, we obtain: 



{e [Y + 3H ( 7 - 1) y - \ Kl iX + ^£ ( 



u 2 _ 4K 

WW 1 WW 2 



X 



X - 2HX - 7 y 



+ 



+ {e [3£~{HY + \hX + £H 2 - ^) X 

4£ 3 H 3 - j*-,) } [>("> + n("> 

e£-fH 2 Y - AK£H 2 + 2H 2 £ 3 h 



+ 



(38) 



When e = 0, Ky' = C^ u ' and this equation reduces to equation (|2l| ) of section 2^ 



It is clear from (^) that KS V > could be conserved in general only if all the five 
coefficients of the equation vanish, a condition that clearly cannot be satisfied by only 
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two arbitrary functions (X and Y). However we are primarily interested in the case of 
perfect fluids, for which F = II = \& = 0. Then for such a fluid either we take adiabatic 
perturbations (£ = 0), or we take it to represent a minimally coupled scalar field [F^ v ' = 
n( y ) = iJ/O) = 0), in which case the entropy perturbation is p£ M = (1 - c^) fiA^ Q. 
We shall now focus our discussion on the latter case, for which equation (|38[) becomes: 



/CM 



y + 3fT ( 7 - 1) y - + ^ - ^) X 



7 



^ + ^ (3c* - 4)] } AC 



+ e\X-jY- 2HX] 



(39) 



Then the condition for K.W to be conserved is that the two coefficients of AM and 
2^ vanish, i.e. only the two arbitrary functions X and y are needed; this justifies our 
choice (|37|) . An equation similar to (|3^) above is readily obtained from (|3^) for adiabatic 
perturbations of a perfect (barotropic) fluid. In a more general case, such as a perfect 
fluid with isocurvature perturbations {£ ^ 0) or an imperfect fluid (as a non-minimally 



coupled scalar field |27|), (37) should be appropriately generalized with the inclusion of 
further arbitrary functions. 

From (|3^) with e = 1, we obtain a pair of first order differential equations for the 
variables X and Y: 



AK \ y 
WW) A 



AH 3 e 2 



ir^r "t- |Tht W<- s 



4) 



and 



X -jY - 2HX = 0. 



(40) 



(41) 



It must be noted that X = X(u, t) and Y = Y(v, t), i.e. both these functions depend on 
the wavenumber v, as well as t: for each v the two equations above determine a pair X, 
Y such that the corresponding harmonic component KP^ of K, will be conserved. These 
two equations can be combined to give a second order equation for X. Substituting for 
7 and in terms of (ft and (ft, we obtain: 



X 



5H + ^ 



X + 



^2 +6H 2 + AHi + 1 (^_ 6K) 



-24# 



7-fQ 



(42) 



We will now apply the above equations to the following two cases commonly used when 
discussing the evolution of energy - density perturbations during an inflationary period: 
a coasting era, characterized by a constant value of the density parameter O and De 
Sitter exponential inflation (see [O, fi^] for a discussion of these models). 



5 Applications 
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5.1 The coasting solution: i{t) = At 



Using the results given in [18], equation ([4^) becomes: 



3 + -tj(v 2 - 5K) 



X = AH [y 2 - 5K) 



(43) 



which we can immediately integrate, giving, as in [12|, three classes of solution: 



a) v 2 < A 2 + 5K 



X(t) = X A t 2+ v + X"t*-W + \AH 



B.2- 



(44) 



where Q = J ^ K A ;i 1 ' + 1 as before and X A , X a are two integration constants. The 
solution for Y then follows from equation (|40|): 



Y(t) = \Q 



6A 2 . 



b) v 2 = A 2 + 5K ; in this case we obtain: 

X(t) = t 2 \X A + X B ln(t)l + 4A 2 t , 

and 

Y(t) = \X B t - QA 2 . 

c) v 2 > A 2 + 5K ; the third case corresponds to damped oscillations: 

X(t) = t 2 [X A cos [Q ln(t)] + X B sin [Q ln(t)] } + 4A 2 t , 



and 



Y(t) = -|Qt{x A sin[Qln(t)] - X s cos [Qln(t)]} - 6^ 2 . 



(45) 

(46) 
(47) 

(48) 
(49) 



Substituting these results into equation ( p7\ ) we can obtain an expression for the conserved 
quantity /C. 

5.2 De Sitter exponential expansion: £(t) = A exput 

This time the second order differential equation for X(t) (|42|) becomes: 



X - 3wX + 



2^ 2 + ^{v 2 - 5K) exp(-2wi) 



X = 4w[v 1 -5K)w 



(50) 



Again, as in the previous example we can integrate this equation to give three classes of 
solution: 
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a) v 2 - 5K < : 



X(t) 



1 



X A exp(2C^/p) + X B exp(-2C7p) 



+ 



C 2 jv 2 -5K) 



(51) 



where C 



\ w <i£i " ) p = exp(— 2ioi) and X" 4 , X s are two integrating constants. The 
corresponding solution for K is: 



y(t) 



^ {x 4 ((1 - p) exp(2C» - ^ exp(-2C v ^) 



+ X B ((1 - p) exp(-2C v ^) + ^ exp(2C») } 



+ 2p-(„2- 51f)(l -p) 



(52) 



b) v 2 — 5K = ; in this case we obtain: 



and 



X(t) = X A exp(-ip) + X B exp(-p) , 



Y(t) = H [x^( p - 2) exp(-ip) + 2X B (p - 1) exp(-p) 



(53) 
(54) 



c) v 2 - 5K > : 



and 



X(t) 



Y(t) 



l 



X A cos(2Cvp) + X u sin(2C^p) 



+ 



C 2 (^ 2 -5A~) 



2cj 



{X^ ((1 - p) cos(2C v ^) + sin(2Cv^) 



+ X B ((1 - p) sin(2CVp) - % cos(2Cyp)) } 



+ 



2C 2 /„2 



(^ 2 -51C)(l-p) 



(55) 



(56) 



Substituting for X(t) and Y(t) into equation ( p7|) we can obtain an expression for the 
corresponding conserved quantity /C. 



6 Transfer functions 

In the previous two sections we have illustrated a method to obtain a conserved quantity 
via two auxiliary functions X(i/,t) and Y(v,t). Without going into further details, 
in the following we shall sketch how to use K,^ v > to connect the amplitude of perturbations 
at different epochs. 

Typically, one is interested in connecting an early universe (inflationary) era with 
the present matter dominated era. For this reason, and for the sake of simplicity, we 
will consider eras in which we either have a perfect fluid with adiabatic perturbations, 
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or a minimally coupled scalar field. In both these cases, the evolution of perturbations 
is given by a system of equations like those of section |2^, i.e. 



AM = fl (t)AM + b(t)zM , (57) 
i:M = C (()AM + d(t)zM , (58) 

where at least one of the functions a, b, c, d depends on v [we remind the reader that 
the exception is the dust case, for which C is conserved at all scales, even in an open 
universe, see fl2l[)1, and all depend on the equation of state valid in the given era. 
We can write equation (p7|) in the form 



£(„) = i)AM + K2 (zy, t)2(") , (59) 

where «x an d K 2 depend on X(t) and Using ([59]) we then obtain an equation of 

the same form as ( p9| ) , and requiring the vanishing of the coefficients of this equation we 
get ([!(]) and (^l]), which can then be integrated (at least in some cases) to give us the 
functions X and Y. With these results in mind, while at the same time using fl59| ) to 
substitute for 2M*> in ([57]), we obtain a first order equation for AM : 

AM =S(t)AW+6(t)2W , (60) 

and since /CM = const, this equation can always be integrated (at least in principle). 
Then, neglecting a decaying mode, this will give 

AM(t) = F{v,t)l& ) , (61) 



or an integral of the form (22) in the general case. If we want to connect two different 
eras with different equation of state, we have to repeat the procedure sketched above for 
the two different cases. Then using one obtains 

A^(t 1 ) = T(u,t 1 ,t 2 )A^(t l ), (62) 

where T(v,ti,t2) = F(y, t2)/F{v, t\) is the required transfer function. 

Thus we have obtained an algorithm through which, provided that the two auxiliary 
functions X and Y can be explicitly determined, the evolution of the density perturbation 
given, for any scale fixed by v, by a first order equation with the conserved quantity 
/CM as source term. Using this, a transfer function is readily obtained, thus generalizing 
the procedure illustrate for example in P, 0] to any v and K = 0, — 1. 



7 Discussion 

In the present paper we discussed under what conditions the curvature variables C and C 
are conserved both in the case of a general imperfect fluid and in detail for a minimally 
coupled scalar field. It was shown that in the general case C is only conserved for 
large scales when the background is taken to be flat or almost flat (i.e. O, ~ 1), while 
for a minimally coupled scalar field an additional condition (^) must be satisfied in 
the almost flat case. We then explored the possibility of defining a gauge - invariant 
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curvature variable fC which would be conserved in a more general set of circumstances 
for all scales of interest. We showed that at least for some models this possibility can 
indeed be realized and illustrated this with the example of a minimally coupled scalar 
field evolving either in a coasting phase (f2 = const) or during a De Sitter inflationary 
era. We also discussed how to use fC to derive a transfer function in order to connect the 
amplitude of perturbations at different epochs. 

Although we restricted our attention to rather special examples, all the equations 
have been derived in general, so this prescription can be applied to a wider variety of 
models both inflationary and non - inflationary. In particular, applications to generalized 
gravity theories will be of great interest in view of the recent growing activity in such 
theories and the generation and evolution of perturbations during a possible inflationary 
phase and their final spectrum at second horizon crossing time deserve special attention. 
These issues will be discussed in a future paper 
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